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CREEPING FLOW SOLUTION OF LEIDENFROST 


BOILING WITH A MOVING SURFACE 

"by Kenneth J. Baumeister 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

and Glen J. Schoessow 
University of Florida 
Gainesville, Florida 

ABSTRACT 

The total vaporization time of water drops in Leidenfrost "boiling 
with supporting plate velocities up to 15 ft /sec has been correlated 
with a standard deviation of 4.5 percent. The analytical model combines 
stagnation and Couette flow profiles. In this paper, the vaporiza- 
tion time for a drop film boiling on a moving surface is related to 
the equivalent vaporization time for a smaller drop placed on a 
stationary surface. 


INTRODUCTION 

The liquid metal space power Rankine System is presently being 
considered as a possible auxiliary power source for space applica- 
tion. In the mist flow regime of its boiler, spiral inserts centri- 
fuge entrained liquid droplets to the heated tube wall to provide 
droplet free vapor to the turbine. At the wall, the liquid droplet 
undergoes film boiling which resembles the classic stationary 
Leidenfrost phenomena, except that in this case a relative velocity 
exists between the drop and the heated tube wall. The present paper 
considers the effect of this relative velocity on the vaporization 
times of the liquid droplets in Leidenfrost film boiling. 

If liquid drops are placed upon a sufficiently hot surface, 
vapor will be generated at the underside of the drop at a rate adequate 
to support the drop (see Fig. 1(a) and (b)). This phenomenon is 
usually referred to as stationary Leidenfrost film boiling. Refer- 
erence 1 considered the Leidenfrost phenomenon on a spinning wheel 
(see Figs. 1(c), (d), and (2)). It showed that the velocity of the 
heated wheel surface had a significant effect on the vaporization 
time of the drop ( see Fig. 3) . 

A model will be developed herein to predict the vaporization 
time of drops on moving surfaces. This model is an extension of the 
authors’ earlier model for Leidenfrost boiling off of stationary sur- 
faces (ref. 2). Table I lists the principal analytical results of 
reference 2 in dimensionless form for a flat disk model as pictured 
beneath the table. 
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Leidenfrost "boiling from stationary surfaces has been under 
extensive study quite recently (refs. 1 to 10). Reference 10, by 
Wachters , contains a comprehensive literature summary of some very 
* early, as well as recent, analytical and experimental work in this 

area. For the present paper; however. Bell's and Hoffman's dis- 
cussion in references 7 and 8, and Bell's educational monograph, 
reference 9, contain pertinent discussions of some of the assump- 
tions and limitations of the stationary Leidenfrost model pre- 
sented in reference 2 upon which the present model is built. 

Recently, Poppendiek et al. (ref. 5) have developed a prac- 
tical correlation for the heat transfer coefficient for helical fog 
flow in a Rarikine boiler using a stationary Leidenfrost model. The 
present analysis will provide a velocity correction factor which 
will account for the relative velocity that exists between the 
drop and the supporting surface. 

LIST OF SYMBOLS 


A 

A* 


B 


1,2 


' 1 , 2,... 8 


C 

F 


P 


n 


f 


g 

Sc 

h s 

hrp 


area of drop 

dimensionless area of lower surface (see table I) 

constants, see table II 

arbitrary constants 

specific heat at constant pressure 

dimensionless function defined by equation (3) 

dimensionless function defined by equation (5) 


radiation factor, 



(reference 1) 
coefficient of gravity 

gravitational constant (conversion constant between 
mass and force units) 

heat transfer coefficient on stationary surface 
total heat transfer coefficient 
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h* _ dimensionless heat transfer coefficient (see table I) 
(conduction mode only) 

h r radiative coefficient, 

h s h s evaluated at V/2 

K constant, equal to l/2(l - 36 2 C Q /u) 

k thermal conductivity 

L length, see figure 4 

l average drop thickness 

Z* dimensionless average drop thickness, see table I 

n correction factor 

P static pressure 

P Q reference pressure 

P„ Prandtl number, 

r k 

r radius of drop 

T temperature 

T w wall temperature 

T g saturation temperature 

AT temperature difference, T w - T s 
t time 

t* dimensionless time, see table I 
U velocity of plate relative to drop 

U* dimensionless velocity defined by equation (4) 
vapor velocity in x direction 
volume of drop 



V 



4 


V* dimensionless volume of drop, see table I 
w vapor velocity in y-direction 

X half side of model, see figure 4 

x distance coordinate 

y distance coordinate 

z distance coordinate 

p p 

(B constant given by equations (A21) and (A22), also p = 

T see bottom table II 

S vapor gap thickness 

liquid emissivity 

X latent heat of vaporization 

/ 7CL £T\-3 

A* modified heat of vaporization, A* = All + — I 

from reference 6 

p viscosity 

Pg vapor density 

PL liquid density 

a surface tension 

a Boltzmann constant 

0 latent to sensible heat ratio , equation (A27) 

t stream function 

stationary stream function 
Couette flow stream function 

METHOD OF ANALYSIS 

The experimental measured vaporization time of a discrete liquid 
drop in Leidenfrost boiling can be estimated by a direct integration of 
an energy balance on the drop: 
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-Xp L || = h T (V)A(V) (1) 

where the total heat transfer coefficient, 1^, and the drop area A 
are dependent on the volume V of liquid that exists at any time t. 
In equation (l), the total heat-transfer rate to the drop, represented 
hy the right side of equation (l), is set equal to the mass evapor- 
ation rate times the latent heat of vaporization. The problem, of 
course, is to relate hij, and A to the properties of the liquid and 
vapor, plate temperature, and the environmental conditions surround- 
ing the drop. 

The mechanism for energy transfer to the drop is assumed herein 
(see appendix A) to he conduction across the vapor film (in creeping 
laminar flow) and radiation to both the top and lower surface. Con- 
duction and diffusive evaporation from the upper surface and sides of 
the drop were neglected. If the vapor concentration above the drop 
is low as in metastable boiling (ref. 11), diffusive evaporation can 
be extremely important, as Bell (ref. 3 and 9) points out and as was 
verified experimentally by Wachters (ref. 10). However, the temper- 
ature differences between the plate and the drop are assumed to be 
very high in this problem. Also, the conduction under the drop is 
much greater in this problem because of velocity effects. Thus, the 
assumption concerning diffusive evaporation will be less restrictive 
in this problem. 

In addition, the shape of the drop as described by the expres- 
sions in table I will be assumed independent of velocity. Photo- 
graphs in reference 1 show that the drop shape depends on velocity. 
However, for the sake of mathematical simplicity, these deviations 
will be neglected. 

With the latter assumptions, only the heat transfer coefficient 
becomes velocity dependent. Thus, to determine the vaporization 
times as a function of plate velocity, the following steps are re- 
quired 

(1) Determine the effect of velocity on the heat transfer 
coefficient. 

(2) Substitute the heat transfer coefficient and the geometric 
relationships from table I into the energy balance, equation ( 1) , 
and integrate to find the vaporization time of the drop. 

The latter step is performed in the next section while the first 
step, that of determining the heat transfer coefficient, is performed 
in the appendix. 

To simplify the derivation of the heat transfer coefficient, the 
drop is represented by an equivalent square shape as shown in fig- 
ure 4. The area of the square is set equal to the circular areas 
listed on table I. A constant gap thickness is assumed beneath the 
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drop. No flow is assumed in the z direction; however, later a cor- 
rection is made for the flow in the z direction (see eq. (A 33 ) ) . 

The physical properties are assumed constant and evaluated at the film 
temperature. Also, the medium surrounding the drop is assumed to be 
saturated vapor and the vapor pulled beneath the drop by the action 
of the moving plate is assumed to have the same temperature profiles 
as that which exists beneath the drop. With these assumptions, the 
heat transfer coefficient derived in the appendix, equation (A34), 
is 


h* - - --y— (2) - (A34) 

1 + F 2 ) 1 ' 2 - £) 

where F as given in equation (2) as 

F = (2K - l)(P r $) 1 /2 u* ( 3 ) . ( A26 ) 

where K is an integration constant to be evaluated by experiment and 
the dimensionless velocity U* given in equation (3) is 


U* 



U 


I72 




\ f° s A 1 / 2 

- 11 g — l * 

) V p c g / J 


(4) - (A28) 


Here, hg represents the dimensionless heat transfer coefficient 
to a stationary drop in Film boiling. The expressions for hg are 
given in table I. The denominator in equation (2) represents a the- 
oretical correction factor which accounts for the effect of velocity 
on the heat transfer coefficient. The theoretical F factor con- 
tains the velocity dependent terms while U* represents the char- 
acteristic dimensionless velocity. 

Empirical Velocity correction 

For large velocities, the theoretical heat transfer coefficient 
is proportional to the half power of the velocity as a combination of 
equations (3) and ( 4 ) will show. This seems to be incorrect, at 
least for the range of data given in reference 1. Consequently, it 
was necessary to modify F in an empirical fashion. A power factor 
n is now introduced into equation (3) to give 

F n = ( 2K - l)(P r $) 1 / 2 U* n 


(5) 
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This correction factor will he carried into the next section' to 
solve for the vaporization time. 

ANALYSIS OF VAPORIZATION TIMES 

Equation (l) can he conveniently written in dimensionless form. 
Expressing hrp, V, A, and t in equation (l) in terms of the dimen- 
sionless forms (listed in table I) gives 

-dV* = h*A* dt* (6) 

The total heat transfer coefficient h^ has been replaced by the 
theoretical h*, equation (2), which considers only the conduction 
mechanism. The relationship between hrf and h* is 

= ^ (7) 

where f represents a radiation factor. The factor f is defined 
in the list of symbols and derived in references (l) and (2). This 
factor considers radiation to both the top and bottom of the drop. 

The f factor was absorbed into the definition of t*. 

Equation (6) can now be integrated for the various domains of 
volume V* defined in table I. Because of the complexity of equa- 
tion (2), a closed form solution could not be obtained. However, by 
evaluating the function F at a reference volume of V*/2 in the 
volume range of interest, a simple closed form expression for the 
vaporization time was obtained. The reference volume technique, 
quite surprisingly, correlated the experimental data better than a 
number of more complicated integration techniques. 

Substituting the value of h* from equations (2) and (5) along 
with h|, 7* and A* from table I into equation ( 6) and performing 
the integration for the various volume domains shown in table I gives 
the results tabulated in table II. 

The equations for the vaporization times, t*, in table II are 
of the same functional form as that in table I for the stationary 
case. In this case, however, V* is replaced by a velocity- corrected 
pseudo-volume V 4- given in the last column in table II. For the 
case when U = 0, the dimensionless velocity U* is zero. Thus, 

F n is zero and the pseudo volume V 4- is identical to V*. 

The integration for V* greater than 0.8 requires that the 
integration be broken up into two ranges since the forms of equa- 
tion (6) are different in each volume regime. Consequently, for 
V* greater than 0.8, the integration is from V* to 0.8 and from 
0.8 to zero. A similar approach is used for V* greater than 155. 
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Breaking up of the integration leads to the constants B^ and Bg 
listed below table II. Of course, a different reference volume is 
used for each volume range. 

SIGNIFICANCE OF PSEUDO VOLUME V+ 

When a drop of liquid with dimensionless volume V* is placed 
on a stationary plate, the dimensionless vaporization time can be 
estimated from the formulas given in table I. However, if the plate 
were given a velocity with respect to the drop, the actual vapori- 
zation time would be much smaller than the vaporization time pre- 
dicted from the stationary equations, as deduced from figure 3. 

The general expressions for stationary vaporization time, how- 
ever, can still be used in the case where a relative velocity exists 
between the drop and the plate. Equations in table II indicate that 
a smaller pseudo dimensionless volume V + if placed on a stationary 
surface would evaporate in the same amount of time as a larger real 
volume V* which is placed on a moving surface. 

COMPARISON OF EXPERIMENT TO THEORY 

The experimental data of reference 1 was correlated by opti- 
mizing the value of K in equation (3) and n in equation (5) 
to give the best fit of the data. The integration constant K was 
1. 3 while the velocity correction parameter n was 0. 4. As seen in 
figure 5, these values correlate the data quite nicely with a 
standard deviation of only 4.5 percent. 

The value of n equal to 0.4 indicates that the heat transfer 
coefficient is a function of velocity to the 0.2 power. This agrees 
closely with the experimentally determined values of 0. 171 and 
0.227 listed in reference 1. 

The major advantage of the present semi-theoretical correlation 
as compared to the purely empirical correlation given in reference 1 
is the possibility that the dimensionless groups and constants would 
be universal and apply to all fluids. Some preliminary data taken 
at the University of Florida by Charles Wood for Ethanol, Carbon 
Tetrachloride and Benezine indicate that the present correlation is 
applicable to a wide range of fluids. 

CONCLUSIONS 

The velocity effects on the vaporization time can be accounted 
for by the introduction of a pseudo volume V + using the property 
groups and constants predicted by a relatively simple flow model. 



APPENDIX - HEAT TRANSFER COEFFIC IEN T 


The heat transfer coefficients from the plate to the drop 
obtained by solving the momentum and energy equations for flow 
heat transfer in the vapor gap beneath the drop. 

are 

and 

For the assumptions noted in the body of the report (see 
Method of Analysis Section) the governing differential equations are 

Momentum : 



0 

oP ( o 2 u d^u\ 

= ' g '^ + V + 

(Al) 

0 

cip ( Aw 

“ - g ° ¥ + ip; &) 

(A2) 

Continuity: 

3u 3v . 

S + 5? = ° 

(A3) 

Energy: 

. . - T s> 

<1 = + 5 

(A4) 

Boundary Conditions : 



y = 0 

u = U w = 0 T = T w 

(A5) 

y = S 

u = 0 w = w(5) 1 = 1,, 

D 

(A6) 


x = +X P = Pq at y = S 

(A7) 


In boundary condition (A5), the drop is fixed and the plate is assumed 
to move with velocity U, as might be approximated by the experiment 
depicted in figure 2. 


Static Force Constraint: 
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Interface Energy Balance (Neglecting Radiation) : 


-p g Aw(5) 


k(T w - T s ) 

5 


(A9) 


The differential energy equation has heen replaced By the simple 
conduction equation (A4) ; however, a correction for the convection 
heating of the vapor will he made later using the results of refer- 
ence 6. Equations (A8) and (A9) are additional and necessary con- 
straints which couple the energy and momentum equations. Equation (A8) 
requires that the integral of pressure Beneath the drop Be equal to 
the drop’s weight. Equation (A9) says that latent heat release at the 
drop surface must Be Balanced By the heat conduction across the gap. 

The Bulk liquid is assumed to Be at the saturation temperature. 

Boundary condition (A7) deserves some special comment. Using a 
uniform gap model to predict the heat transfer coefficient will re- 
quire that the pressure at the front (-X) and the rear (+X) Be dif- 
ferent. As will Be shown later this pressure difference is extremely 
small. Because of the complicated flow pattern around the drop 
(see fig. 2(d)) it is conceptually possible that such a pressure dif- 
ference could occur, although such a pressure difference can Be con- 
sidered pseudo. A model could Be developed in which the pressure is 
equal at the front and rear By using the slider Bearing model (ref. 12, 
page 98) . In this case the gap thickness would Be a linear function 
of position in which the slope could Be determined By some added 
constraint or experiment. However, the uniform gap model was chosen 
Because of its simplicity and Because it does lead to a good corre- 
lation. 


Momentum equations 

The momentum equations are easily solved By use of the stream 
functions: 




(A10) 


w = + 


(All) 


Using these definitions the governing momentum equations can Be com- 
Bined into a single equation with continuity identically satisfied 
(ref. 12, page 59): 


4 3xW 


idx 


T|f = 0 


(A12) 
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The problem now is simply to find a stream function which satis- 
fies equation (A12) and the boundary conditions. To find \|r, assume 

♦ = + 8 (*>y) + t c (y) (Ai3) 

The stream function \|r g is the stationary stream function which would 
’• result if a stationary plate supported the drop, while i|r is a 

stream function representing Couette flow. They are superimposed in 
order to satisfy the boundary conditions. 

Assuming that the variables are separable, gives 

^s = x ( c i + C 2 y + C 3 y 2 + C 4 y 3 ) (A14) 

while >|r can be found by direct integration of equation (A12) 

^c = C 5 + C 6 y + C 7 y2 + C 8 y3 ( A15 ) 

Combining the two stream functions, equations (A14) and (A15) in 
equation (A13), determining the velocities from equations (A10) and 
(All), and applying the boundary conditions in equations (A5) and (A6) 
yields 

u = 3xp 2 (yS - y 2 ) + u|l - ^ - (1 - 2K) (A16) 

w = -p 2 ^| 5y 2 - y 3 ^ (A17) 

Not all the constants C's were determined from the boundary 
conditions. The undetermined constants have been relabeled p 2 and 
K (see monemclature) . The constant p 2 and the gap thickness 5 
will be determined later by applying equations (A8) and (A9). The 
constant K will be discussed later in this appendix. 


Pressure Distribution 


Knowing the velocities under the drop allows us to determine the 
pressure distribution beneath the drop. Substituting the velocity 
distribution back into the momentum equations gives 


= _ e 2 u(56 - 6y) 
by g c 


(A18) 


(A19) 
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Integrating the above equations subject to the boundary conditions 
gives 


P - Pr 


0 


= ^ fix 2 - X 2 ) - (Sy - y 2 )j - : £ (X - X) 

g c L J 


(A20) 


Later, K will be taken as 1.3 (reasons to be given later). This 
will lead to a pressure slightly below atmospheric (approximately 
0.001 atmospheres) at the leading edge of the drop. 

Coupling Momentum and Energy 


Substituting the expression for pressure into the static pres- 


sure balance, integrating, and solving for 

(Pr. 


gives 


e 2 = 


p g) gV + U(4K - 2) 


8(j.X‘ a 


25 2 X 


(A21) 


where L was set equal to 2X, that is, the drop was assumed square. 

The parameter can also be determined by the interface energy 

constraints. Evaluating w at y = 5 in equation (A17), substituting 
this value of w into the interface energy balance, equation (A9), 
and solving for 0^ gives 


2 2k 


(T w - T s ) 


P 4 * = 




(A22) 


The momentum and energy equations are coupled by equating equa- 
tions (A21) and (A22) to give after a little rearranging 


& 4 + 8uX%(2K - 1) & 2 = 16qk £TX 4 

( p L " Pg)g v (Pl - Pg)PgVAg 


Solving for the gap thickness gives 

\l/4 


6 = 


f 16uk Z^T . 

((pl - P g )P^* < f ] 


sr v 


F^ + 1 - F 


1/2 


(A23) 


(A24) 


F = (2K - 1) U 


t 


xfl l/2 

(PL - p g )k ATgvJ 


PPgA* 


(A25) 


or identically after a little manipulation 
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i. 


V' 


F = (2K - IJCP^) 1 / 2 ^ 

where P r is the Prandtl number (see symbols) and 


0 = 


A* 


c p z*r 


(A26) 


(A27) 


u* = 2 

f-HT 


u 


\Pg / \PlS / J 


41 

L. \^g 


(A28) 


(see eq. (A31) for right hand expression) 

The latent heat. A, in equations (A24) and (A27) has been replaced 
by an effective latent heat which is defined in the list of symbols 
and derived in reference 6. This accounts for convection terms in 
the energy equation in an empirical manner. 


Length, Area, Volume Relations 


The length X and the volume V, appear in equations (A24) 
and (A28). These parameters can be related to the geometrical forms 
in table I by the following manipulations 


V = A l = (2 X) 2 Z 


Therefore 




which now can be used in equation (A28). Also, 


V 16V 16V* 



Heat Transfer Coefficient 


(A29) 


(A30) 


(A31) 


The heat transfer coefficient in dimensionless form for pure 
conduction across the gap is 
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k 


k3 ( p L " p g) p L 1 ^ 2p gS 3 ^ 2A * 

1/4 

p. AT0 A / 2 g c 1 /2 





+ 1 



(A32) 


where the value of 5 from equation (A24) was substituted into the 
first part of equation (A32) and V and X in equation (A32) were 
related "by equation (A31). The non-dimensionalizing term used in 
the heat transfer coefficient is the same one listed in table 1 for 
the stationary drop. 

For the case of zero plate velocity, equation (A32) reduces to 


h* = 



(U = 0) 


(A33) 


As can be seen in a comparison to the heat transfer coefficient in 
table I, (V*/A*)V 4 is within twenty percent of the stationary heat 
transfer coefficient which was derived for a symmetric circular drop. 
Thus, the effect of flow in the z direction beneath the drop shown 
in figure 5 is compensated for empirically by replacing (V*/A*)V 4 
by h* which yields 



(A34) 


This correction, however, is exact in the limit of zero plate velocity. 
The denominator represents a velocity correction factor to the sta- 
tionary heat transfer coefficient, h*. 


K-Parameter 


The value of F can not he determined as yet since it is a 
function of the unknown parameter K. 

As seen from equation (A34) , heat transfer coefficient increases 
for increasing values of K. At K = l/2, however, h* equals h* 
and the velocity has no effect on the heat transfer coefficient. 

Thus, K must be chosen such that K > l/2. 
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In addition to affecting heat transfer, the value of K directly 
affects the velocity distribution, as given by the second term in 
equation (A16) . This term is shown graphically in figure 6. Since 
the moving plate is the cause of the real or pseudo pressure differ- 
ence across the drop it is reasonable that no back flow occur; thus, 

K <L 1. Consequently, K appears to be bracketed in the range 
1/2 < K < 1. 

The value of K was determined experimentally by correlating 
the data of reference 1. The experimentally determined K was 1.3, 
as discussed in the body of this report. It is quite remarkable that 
the optimal value of K came so close to the upper theoretical limit, 
as shown in figure 6. 
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Stationary Drops 



Figure 2 . - Cross section of wheel assembly. 
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Figure 6. - Coutte flow profiles as a func- 
tion of K (second term in eq. (A16)). 
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